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Abstract 



We study the behaviour of a radiating star when the interior expanding, shearing fluid 

CT 1 . 

particles are traveling in geodesic motion. We demonstrate that it is possible to obtain 
new classes of exact solutions in terms of elementary functions without assuming a 
separable form for the gravitational potentials or initially fixing the temporal evolution 
of the model unlike earlier treatments. A systematic approach enables us to write the 
junction condition as a Riccati equation which under particular conditions may be 
transformed into a separable equation. New classes of solutions are generated which 
allow for mixed spatial and temporal dependence in the metric functions. We regain 
particular models found previously from our general classes of solutions. 



X 1 Introduction 

In an astrophysical environment, it is likely that a star emits radiation and particles 
in the process of gravitational collapse. In this situation, the heat flow in the interior 
of a star should not be neglected; the interior spacetime of the collapsing radiating 
star should match to the exterior spacetime described by the Vaidya 1 solution. Exact 
models of relativistic radiating stars are important for the investigation of the cosmic 
censorship hypothesis and gravitational collapse^. Santos^ formulated the junction 
conditions for shear-free collapse, matching the interior metric with the exterior Vaidya 
metric at the boundary of the star, which made it possible to generate exact models. 
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This treatment enabled us to investigate physical features such as surface luminosity, 
dynamical stability, relaxation effects, particle production at the surface and temper- 
ature profiles for radiating stars in general relativity. De Oliviera et al— proposed a 
radiating model of an initial interior static configuration leading to slow gravitational 
collapse. It had been shown earlier that the slowest possible collapse arises in the case 
of shear-free fluid interiors^. In a recent treatment Herrera et al 7 proposed a relativistic 
radiating model with a vanishing Weyl tensor, in a first order approximation, without 
solving the junction condition exactly. Then Maharaj and Govender- and Herrera et 
al— solved the relevant junction condition exactly, and generated classes of solutions 
in terms of elementary functions which contain the Friedmann dust solution as special 
case. Later Misthry et al— obtained several other classes of solution by transforming 
the junction condition to the form of an Abel equation of the first kind. These exact 
models have proved to be useful in analysing the relativistic behaviour of a collapsing 
objects in the stellar scenario. 

Another useful approach in studying the effects of dissipation is due to Kolassis et 
al— in which the fluid particles are restricted to travel along geodesies. In the absence 
of heat flow, the interior Friedmann dust solution is regained. This particular exact 
solution formed the basis for many investigations involving physical features such as 
the rate of collapse, surface luminosity and temperature profiles. The physical inves- 
tigations include the analytic model of radiating gravitational collapse in a spherical 
geometry with neutrino flux by Grammenos and Kolassis^, the model describing re- 
alistic astrophysical processes with heat flux by Tomimura and Nunes— , and models 
undergoing gravitational collapse with heat flow which serves as a possible mechanism 
for gamma-ray bursts by Zhe et al—. Herrera et al— investigated geodesic fluid spheres 
in coordinates which are not comoving in the presence of anisotropic pressures. Goven- 
der et al— demonstrated that the temperature in casual thermodynamics for particles 
traveling on geodesies produces higher central values than the Eckart theory. The first 
exact solution, that we are aware of, with nonzero shear was obtained by Naidu et al—, 
considering geodesic motion of fluid particles; later Rajah and Maharaj — obtained two 
classes of nonsingular solutions by solving a Riccati equation. Note that the geodesic 
condition arises in several other astrophysical situations including Euclidean stars, with 
the real and proper radius equal, in the absence of dissipation—. 

In this paper we attempt to extend the previous treatments of Naidu et al— and 
Rajah and Maharaj— through a systematic approach by studying the fundamental 
junction condition. Our intention is to show that the nonlinear boundary condition 
may be analysed mathematically to produce an infinite classes of exact solutions. In 
Section 2, we present the geodesic model governing the description of a radiating star 
using the Einstein field equations together with the junction conditions in the presence 
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of anisotropic pressure. We show that it is possible to transform the junction condition 
into a separable equation by placing restrictions on one of the gravitational potentials. 
Two new classes of solutions are obtained in terms of arbitrary functions of the radial 
coordinate and we regain the models found in the past for particular choices of arbitrary 
functions in Section 3. In Section 4, we discuss some physical aspects of the models 
generated. 



2 The model 

In the context of general relativity, the form for the interior space time of a spherically 
symmetric collapsing star with nonzero shear when the fluid trajectories are geodesies 
is given by the line metric 

ds 2 = -dt 2 + B 2 dr 2 + Y 2 (d6 2 + sin 2 



where B and Y are functions of both the temporal coordinate t and radial coordinate 
r. The fluid four- velocity vector u is given by u a = 5$ which is co moving. For the line 
element (JT]), the four-acceleration u a , the expansion scalar and the magnitude of the 
shear scalar a are given by 

if = 0, (2a) 
B Y 

© = 5+2p (2b) 

° = 3(7-1?)' < 2c > 

respectively, and dots denote the differentiation with respect to t. The energy momen- 
tum tensor for the interior matter distribution is described by 

T ab = (p + p)u a u b + pg ab + ir ab + q a u b + q b u a , (3) 

where p is the isotropic pressure, p is the energy density of the fluid, n ab is the stress 
tensor and q a is the heat flux vector. Anisotropy plays a significant role in gravitational 
collapse and affects the mass, luminosity and stability of relativistic spheres; these 
features have been highlighted in the treatments of Chan— and Herrera and Santos 21 . 
The stress tensor has the form 

T^ab = (p r - Pt) \ n a n b - -h ab ) , (4) 
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where p r is the radial pressure, p t is the tangential pressure and n is a unit radial vector 



given by n a = 45". The isotropic pressure is given by 



V = \{Pr + 2p t ) (5) 



in terms of the radial pressure and the tangential pressure. For the line element ([I]) 
and matter distribution (J3J) the Einstein field equations become 
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where the heat flux q a = (0, q, 0, 0) is radially directed and primes denote the dif- 
ferentiation with respect to r. The system of equations fl6a|) -f )6d|) governs the most 
general situation in describing geodesic matter distributions in a spherically symmetric 
gravitational field. These equations describe the gravitational interaction of a shearing 
matter distribution with heat flux and anisotropic pressure for particles traveling along 
geodesies. From ( I6a|) -( l6d|) . we observe that if the gravitational potentials B(t,r) and 
Y(t,r) are specified then the expressions for the matter variables p,p r ,p t and q follow 
by simple substitution. 

The Vaidya exterior spacetime of radiating star is given by 

2m (v 



ds z = - 1 



R 



dv 2 - 2dvdR + R 2 (d6 2 + sin 2 



(7) 



where m(v) denotes the mass of the fluid as measured by an observer at infinity. The 
line element (JZJ) is utilized to describe incoherent null radiation which flows in the radial 
direction relative to the hypersurface S which denotes the boundary of the star. The 
matching of the interior spacetime ([T|) with the exterior spacetime (jTJ) generate the set 
of junction conditions on the hypersurface S given by 



dt 

Y(Rv,t) 
m(v) s 

(Prh 



Rs dv 

Y -U + Y*- Y - 
2 V B 2 

(qBh. 



dv. 



(8a) 
(8b) 
(8c) 
(8d) 



The nonvanishing of the radial pressure at the boundary S is reflected in equation (I8dj) . 
Equation ( 18dj) is an additional constraint which has to be satisfied together with the 
system of equations ( pa|) - ([6d|) . 
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The junction condition (18d|) in the case of shear-free spacetimes was first derived by 
Santos^, and later it was extended by Glass^ to incorporate spacetimes with nonzero 
shear. On substituting (l6b|) and (|6d|) in ( l8d|) we obtain 

2YY + Y 2 -^— + —YT - 2^-YY' + 1 = (9) 

B 2 B B 2 K ' 

which has to satisfied on S. Equation fl9]) governs the gravitational behaviour of the 
radiating anisotropic star with nonzero shear and no acceleration. As equation (J9]) 
is highly nonlinear, it is difficult to solve without some simplifying assumption. This 
equation comprises two unknown functions B(t,r) and Y(t,r). To generate a solution 
we have to specify one of the functions so that the resulting equation is tractable. 



3 New exact solutions 

For convenience we rewrite equation in the form of the Riccati equation in the 
gravitational potential B as follows 



B 



Y Y 2 1 
F 7 + 2YY' + 2YY' 



V V 
+ ~Y' B ~ 2Y' < 10 ) 



Equation ( flOj) was analyzed by Nogueira and Chan— who obtained approximate so- 
lutions using numerical techniques. To properly describe the physical features of a 
radiating relativistic star exact solutions are necessary, preferably written in terms of 
elementary functions. An exact solution was found by Naidu et al— which was singular 
at the stellar centre. A new class of solutions was established by Rajah and Maharaj—, 
which contain the Naidu et al— models, in which the singularities at the centre were 
shown to be avoidable. 

The Riccati equation f llOp . which has to be satisfied on the stellar boundary E, is 
highly nonlinear and difficult to solve. Rajah and Maharaj 18 obtained solutions in an ad 
hoc fashion by assuming that the gravitational potential Y (t, r) is a separable function 
and specifying the temporal evolution of the model. In this paper we demonstrate that 
it is possible to find new exact solutions systematically without assuming separable 
forms for Y(t, r) and not fixing the temporal evolution of the model a priori. If we 
introduce the transformation 

B = ZY' (11) 



Z = ^ [FZ 2 - 1} , (12) 



then equation ffl0|) becomes 



where we have set 



F = 2YY + Y + 1. 
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We observe that equation (fT2l) becomes a separable equation in Z and t, and therefore 
integrable, if we let F be a constant or a function of r only. In other word, (fT2"j) is 
integrable as long as F is independent of £. We emphasize that in this approach we 
have not made any assumption about the separability of the metric coefficients B(t,r) 
and Y(t, r) or restricted the t— dependence. We demonstrate that this approach leads 
to two new classes of solutions in the following sections. 



3.1 Analytic solution I 

If we set F — 1 then the function Y is given by 

Y(r,t) = [R 1 (r)t + R 2 (r)] 2 / 3 , 



(13) 



where R\(r) and fl 2 (r) are arbitrary functions of r. For this case equation ( Tl2|) becomes 

1 r „, , 4) 



[z 2 ~ 1] • 



2[ J R 1 (r)t + i? 2 (r)] 2 / 3 
On integrating (jHJ) we obtain 

_ 1 + /(r) exp [3(Eit + ggj^/ji] 
~ l-ZWexp^^t + ^Vs/^] ' 

where /(r) is a function of integration. Hence from (fTTj) . ( Tl~3|) and ffl5|) we get 



(15) 



2 [ l + /(r) exp [3(i?it + F 2 ) 1/3 /^i] ' 

3 [ 1 - /(r) exp [3(i2it + R^/Rx] 

Therefore the line element (CQ) takes the particular form 



\R[t + gy 



(16) 



ds 2 
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l + /(r)exp[3(i? 1 t + i? 2 ) 1 / 3 / J R 1 ] 
l-/(r)exp [3( J R 1 t + J R 2 )V3/ jRl ] 

4/3/ 



[i^t + i^] 2 2 
[flit + A 2 ] 2/3 



+ [fl 1 (r)t + i? 2 (r)] 4 / 3 (ci^ + sin 



;i7) 



The line element ( JlTl) is given in terms of arbitrary functions Ri(r), i? 2 (r) and /(r) so 
that it is possible to generate infinite number of solutions for different choices of these 
functions. Observe that if we set 

Rt = 0,R 2 = r 3/2 

then the equivalent of (IT7|) is 

ds 2 = -dt 2 + dr 2 + r 2 (d6 2 + sin 2 #d0 2 ) 

which is the flat Minkowski spacetime. As the curvature then vanishes we require 
Ri 7^ and the t-dependence in Y is maintained. 
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It is interesting to see that for particular forms of the arbitrary functions we regain 
models found previously. If we set 

Ri = R 3/2 , R 2 = aR 3/2 

then the line element ( JTTj) reduces to 

l + /(r)exp [3(t + af/z/R] 



ds 2 = -dt 2 +{t+ a y/ 3 { r" 



1 2 



1 - /(r)exp [3(t + a)V3/ij] 



dr 2 + R 2 (d6 2 + sin 




The line element (fl8l) corresponds to the first category of the Rajah and Maharaj— 
models for an anisotropic radiating star with shear. Furthermore note that if we set 
a = and R = r then (fl8l) reduces to 



rfs 2 = _ rft 2 +t 4/3 



l + /(r)exp [3^/ 3 /r] 
1 - /(r)exp[3t 1 /3/ r ] 



dr z + R\d6 2 + sin' 



(19) 



The metric ( Tl9|) was first found by Naidu et al— in their analysis of pressure anisotropy 
and heat dissipation in a spherically symmetric radiating star undergoing gravitational 
collapse. Note that when 

R 1 = r 3 / 2 ,R 2 = 0, f(r) = 
the line element (fT7|) takes on the simple form 

ds 2 = -dt 2 + t 4 / 3 [dr 2 + r 2 [dQ 2 + sin 2 9d<f) 2 )] . (20) 



The line element (|20|) corresponds to the Friedmann metric when the fluid is in the 
form of dust with vanishing heat flux. 

3.2 Analytic solution II 

If we set F = 1 + R 2 {r) then the function Y is given by 

Y(r,t) = R 1 {r)t + R 2 {r), (21) 
where Ri(r) and R 2 {r) are functions of r only. For this case equation (1121) becomes 

^ [Ri + 1] 



The solution of 



2[R 1 t + R 2 \ 
can be written as 



Z 2 - 



[R\ + 1] 



(22) 



l + ff(r)[i? 1 t + fl 2 ]v / ^+ 1 > 1 
V'flf + l U-^W[i?it + J R 2 ]^+ T/Rl . 
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(23) 



where g(r) is the function of integration. Hence from (jTTj) . fl21]) and (|23|) we get 

1 



B 



l + g{r)[R l b + R 2 }\f*%+ i ! R i 



l-g(r)[R l t + R 2 ]V^+ 1 /^ 



Therefore the line element ([1]) takes the particular form 

1 



[R[t + R' 2 ]. 



(24) 



-dtr + 



[R{ + 1] 



l + g{r)[R 1 t + R 2 }V R l +1 / Rl 
1 - g^Rj + R 2 ]^+ l ' R \ 



[R[t + R' 2 ] 2 dr 



-[ J R 1 (r)t + R 2 {r)] 2 {d6 2 + sin 2 fld0 2 ) 



(25) 



in terms of arbitrary functions Ri(r),R 2 (r) and g(r). Therefore it is again possible 
to generate an infinite number of solutions to (Q. As in §3.1 we cannot set R\ — 
became the Minkowski spacetime is results. 

The new class of solutions given above does contain previously known models. Note 
that when 

Ri = R, R 2 = <xR 



the line element fl25l) reduces to 



ds 2 



-dt 2 + {t+af 



R' 2 



[R 2 + 1] 



l-/i(r)[t + a]v / « 5 +T/R 



dr 2 + R 2 (d6 2 + sin 2 



(26) 

where we have defined the new arbitrary function h(r) = g{r)R^ RZ+1 l R . The line ele- 
ment (I26p corresponds to the second category of the Rajah and Maharaj— anisotropic 
radiating stars with shear. Such solutions are difficult to interpret but could play a 
role in gravitational collapse for strong fields. 



4 Discussion 

The simple form of the exact solutions that have been generated in our treatment 
make it possible to study the physical features of the model such as luminosity, rate 
of collapse, particle production, neutrino flux and temperature profiles. In particular 
explicit forms for the causal temperature can be found utilizing the Maxwell-Catteneo 
heat transport equation 

rh a b q b + q a = ~K {h a b V b T + Tu a ) . (27) 

This has been done for special cases by Naidu et al— and Rajah and Maharaj 
The causal temperature is well behaved in the stellar interior, in particular the causal 
temperature is everywhere greater than the acausal temperature. Other choices of the 
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metric functions in our new classes of solutions (ITT)) and (125)) generate similar behaviour 
and highlight the role of inhomogeneity in dissipative processes. Consequently our new 
general class of exact radiating stars is physically reasonable. 

We have generated a new class of stellar models with shear and expansion in geodesic 
motion. All solutions found previously arise as special cases in our treatment. Previous 
analysis were ad hoc and effectively required that one of the metric functions should 
be a separable function. The resulting Riccati equation could then be solved. In this 
work, we did not assume separability of the metric functions and did not initially 
fix the temporal evolution of the model which is different from earlier treatments of 
this problem. Our more general approach enabled us to solve the Riccati equation 
systematically; we found that the Riccati equation could be transformed to a separable 
equation. Two classes of exact solutions were explicitly demonstrated to the junction 
condition. These were shown to contain the Naidu et al— and Rajah and Maharaj— 
metrics. The fundamental reason that new solutions are possible in that we have 
relaxed the condition of separability in the metric function Y(t,r). 
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